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We derive new features of Bose-Einstein correlations by means of Quantum Field Theory at finite 
temperature, supplemented by operator-field evolution approach. The origin of the dependence of 
the measured two-particle correlation function (as well as the so-called correlation radius) on the 
particle mass has received its correct explanation. The lower bound on the particle emitting source 
size is estimated. 


I. INTRODUCTION 

Over the past few decades, a considerable number of studies have been made on the phenomena of multiparticle 
correlations induced by collisions between elementary particles and heavy ions. It is understood that the studies of 
correlations between produced particles, the effects of coherence and chaoticity, an estimation of particle emitting 
source size play the important roles in high energy physics. Notice that produced particles affect the time evolution 
of the system. What this does mean? The effect of quantum correlations between particles is essentially a statistical 
mechanical phenomenon, and therefore cannot be fully described only within quantum field theory (QFT) in closed 
form. This problem becomes manifest when we define the basic operators in the QFT at finite temperature, and when 
we derive the evolution of them. The common feature is the coexistence of quantum and statistical fluctuations in 
one system composed particles. 

In our previous papers [1,2], we established the master equation in the form of the field-operator evolution (Langevin- 
like [3]) equation, which allows one to gain a better understanding of possible coherent behavior of an emitting source 
of elementary particles. The clear shapes of both two-particle Bose-Einstein and Fermi-Dirac correlation functions 
were observed in the LEP experiments ALEPH [4], DELPHI [5] and OPAL [6] which also indicated a dependence 
of the measured correlation radius on the hadron (tt, K - mesons) mass. Recently, the ZEUS Collaboration at 
HERA reported [7] the results of Bose-Einstein correlations (BEC) between kaons (charged and neutral). The radii 
of production volume for neutral and charged kaons turned out similar. There is no yet the definite explanation of 
the origin of mass dependence above mentioned (see, e.g., [8]). Notice, that the LUND model [9] does not predict 
such a dependence. One of the aims of this paper is to clarify this question and to give an explanation of the origin 
of the dependence of the radius of two-particle source emitter on a particle mass. 

We focused [2,1] on two specific features of the Bose-Einstein correlations clearly visible when the latter are presented 
in the language of the QFT, supplemented by the operator-field evolution. The features we discussed were: 

(j) How the possible coherence of the hadronizing (or deconfined) system (modeled here by the external stationary 

force P that appears in the Langevin-like equations) influences the 2-particle BEC function (72(Q = — p(j)^), in 

which two particles are characterized by their four-momenta and p'^ ? 

{a) What is the true origin of the experimentally observed (J-dependence of the C 2 {Q) in the approach used in [2,1]? 

The only physical meaning of the memory term K(p^), the noise spectral function '(pijPfi), and the so-called coherence 
function a = a{m?', f?, ff) related to them 


P ^ - as e > 0 and 0, (1) 

\uj-K{Pf,)\ n{uj,P) 

have not as yet been investigated carefully in [1]. In formula (1), K{p^) = ipiPu) KPu) with [c(p^), c{p'^)] = 5'^{p^—p'^), 
where p^ = (w,^; n{oj,/3) = {exp [(cu — p)/3] — 1}~^ is the number of Bose-particles in the reservoir characterized 
by the parameter p, (the chemical potential) and the inverse temperature /3. We have already established [1] that 
there is at ^ 0 a finite volume Ho(r) for a two-particle emitter source, and there is a correlation picture. 

The space in which the massive fields settle has its own characteristic length Lch ~ ^ch ■> same sense that the 

massive field with mass nic for the Compton length Ac = rn~^ is served. The appearance of the 5-like distribution 
e) = f —*■ 5(a;) as e —> 0 is a consequence of the finite volume of produced particles (formally, e is given 
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by the physical spectrum). The volume flo{r) is conjugated to a characteristic mass scale Mch, ^o{r) ~ 1/M^f^ in 
such a way that Mch oo when the correlation domain does not exist. Hence, M^h —>■ 0 obeys the condition a ^ 0, 
which means there is no distortion (under the constant force P) acting on the system of produced particles. On the 
other hand, any distortion that disturbs this system for any reason (the coherence function a > 0) leads to a finite 
domain of produced particles {Mch ^ 0). Rather, a strong distortion, a ^ oo, satisfies the point-like region of the 
particle source emitter. 

It is difficult to derive both K{p^i) and using the general properties of QFT. It is worth stressing that 

knowledge of one of these functions leads to an understanding of the other, due of the necessary condition [2] 


J-oa 2 7r UJ-k{Pfj) 


( 2 ) 


which is nothing other than the generalized fluctuation-dissipation theorem. The spectral properties 


\P^{Pp)? 


2 - 2 
— sin(/3u;)|a; — K{pP)\ , 0 < /3 < oo 

OJ 


follow from Eq. (2) with 


I sin(/3u;) 
lim- 

/3^oo TT LU 


S{uj). 


( 3 ) 


II. GREEN’S FUNCTION 

In this paper, we would like to focus on the role of the particle mass, which influences the correlations between 
particles. To solve this problem, one must derive the memory term K{pk) using the general properties of QFT. 

We suppose that we are working with fields that correspond to a thermal field d)(a:) with the standard definition of 
the Fourier transformed propagator F[G{p)] 

F[Gip)] = G{x -y)=Tr {T[$(x)<l>(y)]p4 , (4) 

with pp = /Tre~^^ being the density matrix of a local system in equilibrium at temperature T = under 

the Hamiltonian FI. 

We consider the interaction of $(a;) with the external scalar field given by the potential U. In contrast to an 
electromagnetic field, this potential is a scalar one, but it is not a component of the four-vector. The Lagrangian 
density can be written 


L{x) = d^^*{x)d^^{x) — {nk + C/)$*(x)<I)(a;) 

and the equation of motion is 


+ rn?)^{x) =—J{x), (5) 

where J{x) = U^{x) is the source density operator. A simple model like this allows one to investigate the origin of the 
unstable state of the thermalized equilibrium in a nonhomogeneous external field under the influence of source density 
operator J{x) = U^{x). For example, the source can be considered as (5-like generalized function, J{x) = p p{x, e)<I>(a;), 
in which p{x, e) is a 5-like succession giving the 5-function as e ^ 0 (where p, is some massive parameter). This model 
is useful because the 5-like potential U (x) provides the model conditions for restricting the particle emission domain 
(or the deconfinement region). We suggest the following form: 

J(x) = -S(idk d>(x) -f Jfl(x), 

where the source J(x) decomposes into a regular systematic motion part S(id^) $(x) and the random source Jr(x). 
Thus, the equation of motion (5) becomes 

+vk - E(i(9^)]$(x) = -Jr{x), 
and the propagator satisfies the following equation: 
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[-pI + m^ - E(p^)]G(p^) = 1. (6) 

The random noise is introduced with a random operator r]{x) = —m~‘^ for that the equation of motion looks 

like: 

+m'^[l + r]{x)]}^{x) = -Jr{x). (7) 

We assume that ri{x) varies stochastically with the certain correlation function (CF), e.g., the Gaussian CF 

{r]{x)T]{y)) = Cexpi-z'^ulf^), z = x-y, 

where C is the strength of the noise described by the distribution function exp(—with Lch being the noise 
characteristic scale. Both C and y,ch define the influence of the (Gaussian) noise on the correlations between particles 
that ’’feel” an action of an environment. The solution of Eq. (7) is 

<^{x) = - J dyG{x,y) Jniy), (8) 


where the Green’s function obeys the Eq. 

{V^ +m‘^[l + r){x)]}G(x,y) = S(x - y). 


The final aim might having been to find the solution of Eq. (8), and then average it over random operator y(x). 
Note that the operator M(x) = + ri{x)] in the causal Green’s function 


G{x,y) = 


M{a 


-6(x-y) 


is not definitely positive. However, we shall formulate another approach, where the random force influence are 
introduced on the particle operator level. 

We now introduce the general non-Fock representation in the form of the operator generalized functions 


b{x) = a(x) + -R(x), 


(9) 


b~^ (x) = a~^ (x) + (x), (10) 

where the operators a(x) and a^(x) obey the canonical commutation relations (GGR): 

[a(a;), a(a;')] = [a“''(a;), a^(a;')] =0, 


[a(a;), a^(a;^)] = S(x — x ). 

The operator-generalized functions R{x) and R^{x) in (9) and (10), respectively, include random features describing 
the action of the external forces. 

Both &+ and b obviously define the GGR representation. For each function / from the space 5'(3?oo) of smooth 
decreasing functions, one can establish new operators b{f) and b^{f) 

b{f) = J f{x)b{x)dx = a{f) + J f{x)R{x)dx, 

= J f{x)b+{x)dx = a+{f) + J f{x)R+{x)dx. 

The transition from the operators a(x) and a+(x) to b(x) and b~^(x), obeying those commutation relations as a(x) 
and a“''(a;), leads to linear canonical representations. If both b(f ) and b^{f) create the Fock representation of the 
GGR, one can then find the operator U that obeys the following conditions: 

Uaif)U-^=b{f), 
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Ua+if)U-^=b+if). 

We now move to a simple physical pattern. Let us define the differential evolution (in time) equation, where the sharp 
and chaotically fluctuating function obeying this equation is the main object. Since we are dealing with continuous 
time, we can formulate the stochastic differential equation applied to each (analytical) function under the distortion 
of a random force. In classical mechanics, the stochastic processes in a dynamic system are under the weak action of 
a ’’large” system [10]. ’’Small” and ’’large” systems are understood to mean that the number of the states of freedom 
of the former is less than that of the latter. We do not exclude even the possibility of interplay between these systems. 
When the ’’large” system is in an equilibrium state (e.g., a thermostat state), our method allows one to describe the 
approximation to the distribution of the probability in order to find the physical states in the ’’small” system. On the 
quantum level, the role of such a ’’small” system is played by the restricted region of produced particles, a particle 
source emitter with a definite size, that we study in this paper. 


III. EVOLUTION EQUATION 

Following the idea of classical Brownian motion [11,3] of a particle with a unit mass, a charge g, and velocity v(t) 
in the external, let us say, electric field E, one can write in the following form the formal equation describing the 
evolution (in real time t) of this particle: 


dtv{t) = 7 v{t) + F + gE, 

where F stands for a random force subject to the Gaussian white noise, and 7 is a friction coefficient. 

Referring to [2,1] for details, let us recapitulate here the main points of our approach in the quantum case: the 
collision process produces a number of particles, out of which we select only one (we assume for simplicity that we are 
dealing only with identical bosons) and describe it by stochastic operators b{p,t) and b^{p,t), carrying the features 
of annihilation and creation operators, respectively. The rest of the particles are then assumed to form a kind of heat 
bath, which remains in an equilibrium characterized by a temperature T (one of our parameters). We also allow for 
some external (relative to the above heat bath) influence on our system. The time evolution of such a system is then 
assumed to be given by a Langevin-type equation [2,1] for the new stochastic operator b{p, t) 

idth{p,t) = A{p,t) E F{p,t) E P (11) 

(and a similar conjugate equation for b'^{p, t)). We assume an asymptotic free undistorted operator a{p, t), and that the 
deviation from the asymptotic free state is provided by the random operator R{p, t): a{p, t) b{p, t) = a{p, t)+R{p, t). 
This means, e.g., that the particle density number (a physical number) {n{p,t))ph = + 0(e), where {n(p, t)) 

means the expectation value of a physical state, while {n(p)) denotes that of an asymptotic state. If we ignore the 
deviation from the asymptotic state in equilibrium, we obtain an ideal fluid. One otherwise has to consider the 
dissipation term; this is why we use the Langevin scheme to derive the evolution equation, but only on the quantum 
level. We derive the evolution equation in an integral form that reveals the effects of thermalization. 

Equation (11) is supposed to model all aspects of the hadronization processes (or deconfinement). The combination 
A(p, t) + F(p, t) in the r.h.s of (11) represents the so-called Langevin force and is therefore responsible for the internal 
dynamics of particle emission, as the memory term A causes dissipation and is related to stochastic dissipative forces 
[ 2 ] 


/ +00 

dTK(p,t-T)b(p,T) 

-00 


with K(p, f) being the kernel operator describing the virtual transitions from one (particle) mode to another. At any 
dependence of the field operator K on the time, the function A{p, f) is defined by the behavior of the system at the 
precedent moments. The operator F(p, t) is responsible for the action of a heat bath of absolute temperature T on a 
particle in the heat bath, and under the appropriate circumstances is given by 

/ + 00 7 

■^'if{pp)c(pp)e-^‘^*. 


The heat bath is represented by an ensemble of coupled oscillators, each described by the operator c(p^) such that 
[c(pp), c'^ip'p^)] = <5^(Pp —and is characterized by the noise spectral function f^fpp.) [2,1]. Here, the only statistical 
assumption is that the heat bath is canonically distributed. The oscillators are coupled to a particle, which is in 
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turn acted upon by an outside force. Finally, the constant term P in (11) (representing an external source term in 
the Langevin equation) denotes a possible influence of some external force. This force would result, e.g., in a strong 
ordering of phases leading therefore to a coherence effect. 

The solution of equation (11) is given in S'(3?4) by 

e)]> (12) 

w - K[p^) 

where w in p(u), e) was replaced by new scale top = to/P. It should be stressed that the term containing p{ujp, e) as 
e ^ 0 yields the general solution to Eq. (11). Notice that the distribution p{u}p, e) indicates the continuous character 
of the spectrum, while the arbitrary small quantity e can be defined by the special physical conditions or the physical 
spectra. On the other hand, this p{ujp,e) can be understood as temperature-dependent succession (3), in which 
e —> f3~^. Such a succession yields the restriction on the /3-dependent second term in the solution (12), where at small 
enough T there is a narrow peak at w = 0. 

From the scattering matrix point of view, the solution (12) has the following physical meaning: at a sufficiently 
outgoing past and future, the fields described by the operators are free and, the initial and the final states of 

the dynamic system are thus characterized by constant amplitudes. Both states, (/?(—oo) and (/9(-|-oo), are related to 
one another by an operator S{R) that transforms state i^(—oo) to state i^(-l-oo) while depending on the behaviour of 

RiPn)- 

V3(-l-oo) = (p{R) = S{R)(p{—oo). 

In accordance with this definition, it is natural to identify S{R) as the scattering matrix in the case of arbitrary 
sources that give rise to the intensity of R. 

Based on QFT point of view, relation (9) indicates the appearance of the terms containing nonquantum fields that 
are characterized by the operators R{p^). Hence, there are terms with R in the matrix elements, and these R cannot 
be realized via real particles. The operator function R{pfi) could be considered as the limit on an average value of 
some quantum operator (or even a set of operators) with an intensity that increases to infinity. The later statement 
can be visualized in the following mathematical representation [2]: 

R{Pn) = \/aS(p^,p^), 5(p^,p^) = {d^{p^l)d{p^))f^, 

where a is the coherence function that gives the strength of the average S(p^,p^). 

In principal, interaction with the fields described by R is provided by the virtual particles, the propagation process 
of which is given by the potentials defined by the R operator function. 

The condition Mch 0 (or Ho('c) ~ -> oo) in the representation 

lim S(p^,p') = lim no{r) n{id, (3) exp{-q'^/2) -^n(w,/3), 

with [2] 

^o{r) = \rorzr‘l 

means that the role of the arbitrary source characterized by the operator function R(pfi) in b(p^) = d(p^) -I- R(p^) 
disappears. 

IV. GREEN’S FUNCTION AND KERNEL OPERATOR 

Let us go to the thermal field operator <I>(a;) by means of the linear combination of the frequency parts (jy^{x) and 
<j)-{x) 

^(a;) = [(j)+{x)+(j)-(x)\ (13) 
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composed of the operators and b'^iP/j.) as the solutions of equation (11) and conjugate to it, respectively: 


(j) {x) = 


(Pp 


{2ttY2{P + 


h+{pYe^^y 


One can easily find two equations of motion for the Fourier transformed operators and in S'(5ft4) 


[uj - K{p^)]b{pY = F{p^) + p{u}p, e), (14) 

[u - K~''{pP]b^{pP = F^{pP + p{ujp,e), (15) 

which are transformed into new equations for the frequency parts 4>'^{x) and (j)~ix) of the field operator <i>(a;) (13) 

ido(j3~^(x) - [ K{x - y)<j3~^{y)dy = F{x) + PdoD{x)e~''\ (16) 

JU4 

—ido(j)~{x)— ( K~'~{x — y) (j)~ {y)dy = F~'~{x) + P doD{x) e~'^*. (17) 

Jui 


Here, the field components 4>'^{x) and (f>~ (x) are nonlocalized under the effect of the invariant formfactors K{x—y) and 
K'^{x — y), respectively. In general, these formfactors can admit the description of locality for nonlocal interactions. 
The function D{x) in Eqs. (16) and (17) obeys the commutation relation 

[$(a;),$(y)]_ = -iD{x) 


and can be written [12] 


Dix) = ^e(a;°) 


S(x^) - O(x^) Ji 

o /„2 V V 


where e(x°) and 0(x^) are the standard unit and the step functions, respectively, while Ji(a;) is the Bessel function. 
On the mass-shell, D{x) becomes 


D{x) 





At this stage, it must be stressed that we have new generalized evolution Eqs. (16) and (17), which retain the general 
features of the propagating and interacting of the quantum fields with mass m that are in the heat bath (reservoir) 
and are chaotically distorted by other fields. Eor further analysis, let us rewrite the system of Eqs. (16) and (17) in 
the following form: 

idQ(lP{x) - K{x)-k4)+{x) = f{x), (18) 


(x) - A:+(x) (x)=/+(a;), (19) 

where A{x) -k B{x) is the convoluted function of the generalized functions A{x) and B{x), and 

f{x) = F{x) + P dQD{x)e~^K 

Applying the direct Fourier transformation to both sides of Eqs. (18) and (19) with the following properties of the 
Eourier transformation 


F[K{x)-k(j)'^{x)\ = F[K{x)]F[(j)'^ {x)], 
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we get two equations 


[-/ - K(p^)]^+(p^) = F[f{x)], 


( 20 ) 


[p° - K+{Pf,)]^ (p^,) = F[f+{x)]. (21) 

Multiplying Eqs. (20) and (21) by — K^{p^) and —— K{p^), respectively, we find 

[-p° - k{p^)][p° - k+{p^mp^) = r(p^), (22) 

where 

r(p^) = [p° - k+{p^)]F[f{x)\ - [p° + k{p^)]F[f+{x)]. 

We are now at the stage of the main strategy: we have to identify the field <I)(a;) introduced in Eq. (4) and the field 
<I)(a;) (13) built up of the fields (j)'^ and as the solutions of generalized Eqs. (16) and (17). The next step is our 
requirement that Green’s function G(p^) in Eq. (6) and the function r(p^,/3), that satisfies Eq.(22) 


[/ + k{p^)][-p° + k+{p^)]f{p^) = 1 , 
must be equal to each other, where [12] 


(23) 


G{pk) ^ G(p ,3 ,m ) ~ 


2 2 2 n _ 1 


w? — p^ — ie 


with g being the scalar coupling constant and the one-loop correction of the scalar field ^ << 1/m^ at 1/4 < 
(rn?/p^) < 1. This means we define the operator kernel K(p^) in (14) from the condition of the nonlocal coincidence 
of the Green’s function G{p^) in Eq. (6), and the thermodynamic function r(p^,/3) from (23) in 5'(3?4) 


lim F[G(p^) -r(p^,, 

x—x' r-^Oyr) 


We can easily derive the kernel operator K{p^) in the form 

k{pk} = {vk Fp^k 


= 0 . 


l + g i(p) ,m) [I- 


CO 


(24) 


where 


and 


= 


1 


2 TT 


967r2m2 VVS 


— 1 ) , p ckL mr 


^{p ,m ) = 


1 


96 7rm2 


i\ 1- 


4m^ TT 


V3, 


, p ~ 4to"^. 


The ultraviolet behaviour at jp^j >> leads to 




-1 


32 k p^ 
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V. CORRELATION FUNCTION 


Out of many details (for which we refer to [2,1]) the 2-particle BEC function for identical particles is defined as 

f(r) ^ 

cm = x{N) ■ ~ Y’, A{rfQ) x{N) ■ [l + (1 + • Q)'), (25) 

J (Pfl) ' J \P^) 

where /(p^,p(.) = ( 6 +(p^) 6 +(pP 6 (p^) 6 (p),)) and /(p^) = (5+(p^)5(p^)) are the corresponding thermal statistical 
averages with b{pfj,) being the corresponding Fourier transformed solution (12). The multiplicity N depending factor 
is equal to x{^) = {N{N — 1))/(A^)^. We have introduced in (25) the factor A{rfQ) ~ 1 -|- X]i>i(^/ ' *5)^ which is 
nothing other than the consequence of the Bogolyubov’s principle of weakening of correlations atlarge distances that 
is characterized by the parameter of weakening of correlations rf relevant to the particle pair correlations. In the 
limit of infinitely small scale rf the function (72 (Q) coincides with that obtained in [2]. The latter statement means 
that one returns to the ideal correlation pattern. 

As was shown in [ 2 ], (note that operators R{p^i) by definition commute with themselves and with any other operator 
considered here), 

/(Pa*>Pm) = + (a+(PA«)a(p)x))(a+(p)x)a(PA*)) + ip'^iPr) + (p'^Pr))^'^ (Pr)R(p'^)^ 

= {C^iP^Pr)) + \R{Pr)?- (26) 

This defines D{p^,p'^) = /(p^,p(j)/[/(P/i) • /(p)j)] — 1 in Eq. (25) in terms of the operators a(p^) and R{Pfi) which 
in our case are equal to 


KPr) 


Hpj^) 

LO - K{Pf,) 


and 


HPf^) 


Pjt^p, e) 

w - KiPf,) 


This means, therefore, that the correlation function C 2 {Q), as defined by Eq. (25), is essentially given in terms of 
p{ujp,e) and the following two thermal averages for the thermostat operators F{p,t) (for details see [13]): 


{F+{p,t)F{j/,t')) = 6^{p-]^) J ^ |V’(p^)|%(w,/3)e+“(‘ 
{F{p,t)F+{]^,t')) = 6^{p- J ^ |V'(Pm)I^ [ 1 -bn(w,/3)]e"“(‘"*'). 


(27) 

(28) 


Notice that with only delta functions present in Eqs.(27) and (28) we would have a situation in which the hadronizing 
(or deconfined) system would be described by some kind of colored noise only because of the presence of n{uj, /3) that 
carries the quantum properties . The integrals multiplying these delta functions and depending on (a) the momentum 
characteristic of a heat bath ^piPn) (^) assumed bosonic statistics of produced secondaries resulting in factors 
n{uj,(3) and 1 + n{ijj,j3), respectively, bring the description of the system considered here closer to reality. 

We now see easily that the existence of BEC (i.e., that C 2 {Q) > 1) is strictly connected with nonzero values of 
the thermal averages Eqs. (27) and (28). However, in the form presented there, they differ from zero only at one 
point, namely for (5 = 0 (i.e., for p^ = p'^). Actually, this is the price we pay for the QFT assumptions tacitly made 
here, namely for the infinite spatial extension and for the uniformity of our reservoir. However, we know from the 
experiments in, e.g., [14,4-7] that (7*2 ((5) reaches its maximum at Q = 0 and falls towards its asymptotic value of 
(72 = 1 at large of Q (actually at (3 ~ 1 GeV/c). To reproduce the same behaviour by means of our approach, we 
must replace the delta functions in Eqs. (27) and (28) by functions with supports larger than those limited to one 
point only. This means that these functions should not be infinite at (3^ = Pn ~ P'n “6 but remain more or less 
sharply peaked at this point, otherwise remaining finite and falling to zero at small but finite values of \Qp\ (actually 
identical to those at which (72 (( 3 ) reaches unity) 


KPu-p'u) ^2o-exp[-(p^-p;,)L^'^(r)(p^-p;,)]. (29) 

Here we replace the (5-function with the smearing (smooth) dimensionless generalized function f2(g = Qr) = exp[—(py^ — 
p'^)L^^''{r){pi, —p^)], where is the structure tensor of the space-time size and it defines the spherically-like domain 
of emitted (or produced) particles. 

We therefore tacitly introduce a new parameter, a 4-vector that has the dimension of length. This defines the 
region of nonvanishing particle density with the space-time extension of the particle emission source. Expression (29) 
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must be understood in the sense that Q{Qr) is a function that, in the limit of r ^ oo, strictly becomes a 5 - function. 
With such a replacement, we now have 


where 




(1 + q;)(1 + a') 


\/^) + 



n(q) = j ■ [}(q), 7 


n{uj,j3)n{uj',(3)' 


(30) 


(31) 


The coherence function a is another very important one that summarizes our knowledge of other than space-time 
characteristics of the particle emission source. Notice that a > 0 only when P ^ 0. For a = 0, we actually find 


1 < C 2 {Q) < x(lV)[l -k 7fl((3r)](l + rf - QP ...), 


i.e., it is contained between the limits corresponding to very large (lower limit) and very small (upper limit) values of 
P. Because of this, a plays the role of a coherence parameter. Ignoring the energy-momentum dependence of a, and 
assuming that o! = a, we get the expression 


C 2 iQ) ^ X{N) ^1 + Xnewim,f3)e I + Xcorrim, (3) e'^'^ 


} (1 + ^/ ■ Q + ■•■)> 


(32) 


where the new intercept function becomes as Xnew = 7 ( 1 ^)/?)/(! + cr )^7 and the new coherence correction in the 
brackets of Eq. (32) carries an additional intercept function Xcorr = 2 ajP)- In fact, since a ^ ex' (because 
w ^ to' and, therefore, the number of states identified here with the number of particles with given energy n(a;) is 
also different), we must use the general form Eq. (25) for C 2 with details given by Eqs. (30) and (31) and with a (1) 
depending on the particle mass and such characteristics of the emission process as the temperature T and chemical 
potential p occurring in the definition of n{uj). Note that Eq. (32) differs from the usual empirical parameterization 
ofC 2 (Q) [15,16,4-7], 

C 2 (Q) ~ (1 + Ae-«')(1 +ag+ ...), (33) 


which is nothing other but the Goldhaber parameterization [17] at a = 0 with 0 < A < 1 being a free parameter 
adjusting the observed value of ( 72 (Q = 0), customary called a ’’coherence strength factor” or ’’chaoticity” of A = 0 
for fully coherent sources and A = 1 for fully incoherent sources; a is a c-number, with Q{Qr) usually represented as 
Gaussian. Returning to the particle mass dependence of C 2 {Q) (the correlation radius, in particular), we find that 
this dependence comes from a-coherence function ( 1 ) containing the operator kernel iF(p^), defined correctly up to 
the second-order of the scalar coupling constant in Eq. (24) within the framework of the QFT. 

The a - representation in (1) must be clarified: In fact, can be broken down into two parts: ^ ‘ 

where is the small massive scale characterized by the time-like scale r°, while Mch is the characteristic mass 
associated with the spatial inverse components rz, rt, ~ (rzr^)~^. Taking into account the properties of the 

distribution [p{lop, e)]^, we can suggest the following replacement e)]^ ^ [^ch p{^y ^)][^ch Pi^Py s)]) where 

the first multiplier is of the order of 0 ( 1 ), while the second reflects the massive scale Pch of the particles production 
region, i.e., ^IhPit^Py^) ~ Pch < m. Thus, 


a 


__ 

_\u - K{pf_,)\^ niuj,f3)_ 


(34) 


Let us return to the problem of the Q-dependence of the BEG. One more remark is in order here: The problem 
with the S{p^ ~p'^) function encountered in two particle distributions does not exist in the single particle distributions 
that are in our case given by Eq. (26), and which can be written as 


hPr) = (1 + a) ■-{PfiyPr)y 

where 5(p^,p^) is the one-particle distribution function for the ’’free” (undistorted) operator a{p^), namely 


( 35 ) 


iPnyPh) = (a^(PM)d(P/x)) = ^0 


•-Kipfz) 


i(w,/3). 
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Notice that the actual shape of f{p^j) is dictated by both n(w) (calculated for fixed temperature T and chemical 
potential p at energy w, as given by the Fourier transform of field operator K (24) and the shape of the reservoir in 
the momentum space provided by and by the 6 - like distribution of external force p{uip,e). On the other 

hand, it is clear from Eq. (35) that {N) = {Nch) + {Ncoh), where {Nch) and (Ncoh) denote multiplicities of particles 
produced chaotically and coherently, respectively. 

The m - dependence of C 2 {Q,m) is essentially given by 


a{rn?, (3) 


n{ui,(3) 


(36) 


which is nothing other but the effective number of Bose-particles in the plane phase-space with the size Lch having 
the mean mass to(w, j3) = mn{ui, /3). Obviously, a ^ 0 as n(u}, j3) ^ oo (Goldhaber parameterization) and a ^ oo as 
n(u;, /3) ^ 0 (that is trivial result with C 2 {Q) — !)■ Notice, that the condition a(m, (3) <l always yields at sufficiently 
strong temperature T, and its validation would be very weak when T ^ 0. In the latter case the heat bath is simply 
absent. 

On the other hand, neglecting the energy-momentum dependence of a {a'), the latter can be estimated within the 
formula 


a ~ 


2-C2(0) +V2-C2(0) 
C2(0) - 1 


where (72(0) = (72(0)/x(N(), and we assume that the BEG function at (5 = 0, (72(0), and the mean multiplicity 
{N) are known from the experiment. Within the formula (36) and the Eq. (32), it is evidently that an increasing 
of the particle mass m leads to an enhancement of the intercept function Xnewijn, P)- Therefore in the case of the 
investigation of the correlations between heavy particles the chaotic coherence effects become negligible (a ^ 0). 

Finally, we present the lower bound on the particle emitting source size: 

1 

1 ] ^ 

7(n)/A(7^“"^ - ij ’ 


ch ^ Lch ^ 


m 


\/n(0) 


where A(F^“"= = - 1, (F^“"= = (G^“"=/x(A), and is the maximal value of the G 2 -function in the vicinity 

of (5 = 0. 

The lower bound on the particle emitting source size carries the dependence of: 

- the particle mass m given by a and a' which are defined within (34) with the kernel operator K{pfi) (24) calculated 
in the QFT; 

- the mean multiplicity factor x(-^) defined within formula (25); 

- the maximal value of the BEG function C 2 ', 

- the absolute temperature of a heat bath and the chemical potential (the 7 (n)-factor dehned in Eq. (31)). 


VI. CONCLUSION 

To summarize: using the QFT, supplemented by Langevin-like evolution Eqs. (16) and (17) to describe hadroniza- 
tion (or deconhnement) processes, we have derived the two-particle BEG function in a form explicitly showing the 
origin of both the so-called coherence (and how it influences the structure of the BEG) and the (^-dependence of the 
BEG represented by the correlation function (72 ((5)- The dynamic source of coherence is identified in our case with 
the existence of a constant external term P in the evolution equation. Therefore, for P —>■ 00 , we have all phases 
aligned in the same way, and (72 (C?) = 1. This is because coherence has already been introduced on the level of a 
particle production source as a property of the fields or operators describing produced particles. It is therefore up 
to the experimenter to decide which proposition is followed by nature: the simpler formula (33) or the rather more 
complicated Eq. (25) in combination with Eq. (30). 

With our approach, it is also clear that the form of C 2 reflects the distributions of the space-time separation between 
the two observed particles. 

Finally, we would like to stress that our discussion is so far limited to only a single type of secondaries produced. 
This is enough to attain our general goals, i.e., to explain the possible dynamic origin of coherence in BEG, the origin 
of the specific shape of the correlation (72 (C?) functions, and to explain the dependence of the correlation radius on 
the particle mass due to coherence function a, as seen from the the QFT perspective. Actually, r decreases with the 
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mass m. It is then plausible that, in describing the general BEC effect, they should be combined somehow, especially 
if the experimental data indicate such a need. 

The final note concerns to that the ZEUS at HERA resulted that r values for pions and kaons are not so different and 
the effects comes from heavier particles. Therefore more precise measurements from different processes are necessary. 
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